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METHOD AND APPARATUS FOR DETERMINATION OF AN 
ACOUSTIC RECEIVER'S POSITION 



1. FIELD OF THE INVENTION 

The present invention pertains to determining the position of an acoustic receiver and, 
more particularly, to determining the position in an apparatus for seismic surveying. 

2. DESCRIPTION OF THE RELATED ART 

Seismic exploration is conducted on both land and in water. In both environments, 
exploration involves surveying subterranean geological formations for hydrocarbon deposits. 
A survey typically involves deploying acoustic source(s) and acoustic sensors at 
predetermined locations. The source(s) imparts acoustic waves into the geological 
formations. Features of the geological fonnation reflect the acoustic waves to the sensors. 
The sensors receive the reflected waves, which are the processed to generate seismic data. 
Analysis of the seismic data tiien indicates probable locations of the hydrocarbon deposits. 

Accurate knowledge of source and sensor positions is important to the accuracy of the 
analysis. In land surveys, accurate positioning is not particularly difficult because 
environmental conditions are usually relatively stable. Sources and sensors can be readily 
positioned where desired and, once placed, they usually do not shift to any great degree- 
Marine surveys, however, are different altogether. 

Marine surveys come in at least two types. In a first, an array of streamers and sources 
is towed behind a survey vessel. In a second type, an array of seismic cables, each of which 
includes multiple sensors, is laid on the ocean floor, or sea bottom, and a source is towed ficom 
a survey vessel. In both cases, many factors complicate determining the position of the 
sensors, including wind, currents, water depth, and inaccessibility. 

In the second type of marine survey, where the array of seismic cables is laid on the 
sea floor, much attention is paid to the positioning of the seismic cables as they are laid. One 
important consideration is the shape of the seismic cables as they are deployed. The catenary 
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shape, or the shape of the seismic cable in the water during deployment, of a seismic cable 
must be known, or projected, if it is to be controlled effectively during deployment. Control 
is needed to optimize the deployment speed and accuracy. Control is also desired to avoid 
tangling the seisnfuc cable with other obstructions, such as other cables or sub-sea devices. 
Remedial action can be taken to avoid such problems and improve the safety of sub-sea 
operations. 

Current techniques apply various modeling techniques to project the shape and/or 
position of the seisuMC cable during deployment. These models consider the physical 
characteristics of the seismic cable (e.g., weight, diameter, etc.) and account for the effect of 
predicted sea currents on tfie seismic cable as it descends to the sea floor. However, such 
methods provide only a model, or projection, of the seismic cable's shape and are predicated 
on a limited knowledge of the sea's properties. 

The present invention is directed to resolving, or at least reducing, one or all of the 
problems mentioned above. 

SUMMARY OF TPF. TNTyRNTf OM 

The invention comprises, in its various embodiments and aspects, a method and 
apparatus for determining a position of an acoustic receiver. Hie method includes 
determining a plurality of acoustic ranges from at least a first signal source position and a 
second signal source position, respectively, to the acoustic receiver, ascertaining a non- 
acoustic constraint on the acoustic receiver's position; and determining the acoustic receiver's 
position from the first and second acoustic ranges and the non-acoustic constraint. The 
apparatus includes at least one acoustic source; an acoustic receiver, and a computing system. 
The acoustic receiver is capable of receiving a plurality of acoustic signals' transmitted by the 
at least one acoustic source from at least two signal source positions. Hie computing system 
is programmed to determine a position of the acoustic receiver from the acoustic ranges 
between the at least two signal source positions and the acoustic receiver and a non-acoustic 
constraint. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention may be understood by reference to the following description taken in 
conjunction with the accompanying drawings, in which like reference numerals identify like 
elements, and in which: 

FIG. 1 depicts the deployment of a subsea seismic cable in accordance with one 
particular embodiment of the present invention showing a representation of a cable catenary 
measurement using a combination of acoustic and angular measurements; 

FIG. 2 illustrates, in a partially cut away view, one sensor module of the seismic cable 
in FIG. 1; 

FIG. 3 represents a method practiced in accordance with the present invention to 
determine the position of an acoustic receiver; 

FIG. 4 and FIG. 5 conceptually illustrate the analytical determination of the position 
of a sensor module as deployed in FIG. 1 in accordance with a first embodiment of the 
present invention; 

FIG. 6 represents one particular implementation of the method in FIG. 3 wherein the 
positions of points on the cable are detCTnined analytically as illustrated in FIG. 4 and FIG. 
5; 

FIG. 7 represents a one particular implementation of the metiiod in FIG. 3 alternative 
to that of FIG. 4 - FIG. 6 wherein the position is determined using a model-based solution; 

FIG. 8 is a block diagram depicting a computing system such as may be used in 
implementing certain aspects of the present invention; 

FIG. 9 re^niesents a metiiod by which the present invention may for dynamically 
determining the shape of a body in accordance with a first implementation of one application 
of the present invention; 

FIG. 10 represents a method by which the shape of a cable may be determined 
dynamically practiced in accordance with a second implementation of the application; 

FIG. 11 represents a method wherein the positions of points on the cable are 
detmnined through application of an inversion algorithm; 

FIG. 12 projects tiie cable position onto the x-y plane (top) and x-z plane (bottom) for 
a single acoustic source experiment simulating an application of the invention to determine 
the shape of a seismic cable as it is deployed; 
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FIG. 13 graphs the cable position error for the single acoustic source experiment of 
FIG. 12; 

FIG. 14 projects the cable position onto the x-y plane (top) and x-z plane (bottom) for 
a three acoustic source experiment simulating an application of the invention to determine the 
shape of a seismic cable as it is deployed; 

FIG. 15 graphs the cable position error for the three acoustic source experiment in 

FIG. 14; 

FIG. 16 projects the cable position onto the x-y plane (top) and x-z plane (bottom) for . 
a three source and dip angle experiment simulating an application of the invention to 
detemiine the shape of a seismic cable as it is deployed; and 

FIG. 17 graphs the cable position error for the three source experiment and dip 
expmmCTt in FIG. 16. 

While the invention is susceptible to various modifications and alternative forms, the 
drawings illustrate specific embodiments herein described in detail by way of example. It 
should be understood, however, that flie description herein of specific embodiments is not 
intended to limit the invention to the particular forms disclosed, but on the contrary, the 
intention is to cover all modifications, equivalents, and alternatives as defined by the 
appended claims. 

nRTAn.Rn nEscRipnoN of the invention 

Illustrative embodiments of the invention are described below. Jn the interest of 
clarity, not all features of an actual implementation are described in this specification. It wUl 
of course be appreciated that in the development of any such actual embodiment, numerous 
implementation-specific decisions must be made to achieve the developers' specific goals, 
such as compliance with system-related and business-related constraints, which will vary 
fiom one implementation to another. Moreover, it wiH be s^preciated that such a 
development effort, even if complex and time-consuming, would be a routine undertaking for 
those of ordinary skill in the art having the benefit of this disclosure. 

FIG. 1 conceptually illustrates a seismic cable 100 being deployed fix>m a vessel 105 
that is also towing a first source 110 and a second source 115. Hie seismic cable 100 
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comprises a plurality of sensor modules 120 - 126 on a cable 130 and temoinating in a 
steerable anchor 145. The sensor modules 120 - 126 are, in the illustrated embodiment, 
evenly spaced apart a known, predetermined distance, along the seismic cable 100. As is 
ihown in FIG. 2 for a sensor module 200, each sensor module includes an acoustic receiver 
205 and contains, in the illustrated embodiment, an angular orientation measurement device 
210 in a housing 215. The acoustic receivers 205 and angular orientation measurement 
devices 210 receive power over leads 220 through the cable 130. The sensor modules 120 - 
126 may also include additional acoustic receivers (not shown) for use in a seismic survey 
that may be either hydrophones (not shown) or geophones (also not shown). 

Note that, in alternative embodiments, the sensor modules 120 — 126 may be spaced 
apart at irregular distances or clustered on the seismic cable 100. Note also that, in alternative 
embodiments, the' acoustic and angular orientation measurement devices may be housed 
separately. TTie angular orientation measurement devices 210 may be, e.g., inclinometers, but 
other types may be used. The sensors 120 - 126 may also house other kinds of sensors, e.g., 
magnetometers or compasses for measuring heading. The steerable anchor 132 may be 
controlled from the surface 150 of the sea 155 (shown in FIG- 1) in conventional fashion as 
the seismic cable 100 descends to the sea floor 160 to help position the seismic cable 100 as 
desired. 

The scenario depicted in FIG. 1 is illustrative only. For example, the acoustic sources 
110, 120 may be towed from vessels other than the vessel 105 from which the seismic cable 
100 is deployed. The acoustic sources 110, 115 may also be deployed from separate vessels 
in some embodiments or from buoys. Similarly, the seismic cable 100 is shown with seven 
sensor modules 120 - 126, even though the invention is not limited by the number of such 
sensor modules. Indeed, the invention may be applied to embodiments in which the seismic 
cable, e.g., the seismic cable 100, includes any number of sensor modules, e.g., the sensor 
modules 120 - 126. Those skilled in the art having the benefit of this disclosure will 
nevertheless appreciate that certain implementation-specific, practical considerations will 
circumscribe the range of the number of sensor modules. 
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As the seismic cable 100 is deployed, enviroiunental conditions, such as wind and 
current, impart forces on the seismic cable 100 and the vessel 105. These forces distort the 
path of the seismic cable 100 along all three coordinate axes jc, y, and z, causing deviations in 
the path. These deviations, in turn, affect the position of the sensor modules 120 - 126 — and, 
hence, the position of any sensors and/or receivers they house — on the sea floor 160. The 
state of the art provides numerous techniques by which these deviations can be projected or 
predicted, but none by which they can be empirically determined. Thus, the state of the art 
fails to provide a technique by which the position of the sensor modules 120 - 126 can be 
known. The present invention, however, remedies this deficiency. 

FIG. 3 represents a method 300 for determining a position of an acoustic receiver, 
e,g., an acoustic receiver 305, in accordance with the present invention. The method begins 
by first determining (at 310) a plurality of acoustic ranges firom at least a first signal source 
position (e.g., the position of the first acoustic source 110) and a second signal source position 
(^.g., the position of the second acoustic source 115), respectively, to the acoustic recover 
305, A non-acoustic constraint (e.g., the angular orientation of the sensor module 121) on the 
acoustic receiver's position is also ascertained (at 320). The acoustic receiver's position is 
then determined (at 330) from the first and second acoustic ranges and the non-acoustic 
constraint. 

The invention admits variation in the practice of the method 300. For instance, the 
non-acoustic constraint can be ascertained (at 320) prior to determining the acoustic ranges (at 
310). A wide variety of non-acoustic constraints may also be employed, hi the illustrated 
embodiment, the non-acoustic constraint is the respective angular orientation of the sensor 
modules 120 * 126 as measured by the respective angular orientation measurement devices 
210. However, as will become apparent from the discussion below, other non-acoustic 
constraints may also be employed. The determination of the acoustic receiver's position may 
also be implemented in various ways. Two such implementations are discussed in further 
detail below. The present invention may also find many applications. The illustrated 
embodiment employs the invention not only to determine the positions of the acoustic 
recovers, but to also determine the shape of the seismic cable 100. In one implementation. 
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the shape of the seismic cable 100 is detennined dynamically as the seismic cable 100 is 
deployed. 

The deteraiination of the acoustic ranges (at 320, in FIG. 3) can be implemented in a 
variety of ways. Two alternative embodiments are discussed further below. The first is an 
analytical-based solution illustrated in JIG. 4 - FIG. 6. This analytical determination 
employs the analysis represented by Eq. (1) - Eq. (11) and the discussion that follows. The 
second embodiment is a model-based approach illustrated in FIG. 7 and employing die 
analysis lepiesented by Eq. (1) - Eq. (3) and Eq. (12) - Eq. (22). This approach applies an 
inversion algorithm to correct positions initially determined by a conventional model. Tlie 
inversion algorithm iteratively applies a non-linear inversion represent by Eq. (22) to the 
predicted positions using the dynamic angular orientations. However, the invention admits 
wide variation in performing the determination (at 620). Additional alternative embodiments 
may become apparent to those in the art having tiie benefit of this disclosure. 

Referring again to FIG. 1, as the seismic cable 100 is deployed, the acoustic sources 
110, 115 emit acoustic signals 165, Note that, in a seismic cable 100 including N sensor 
modules 120 - 126, there are 3 x iV degrees of freedom. The seismic cable 100 enters the sea 
155, in the illustrated embodiment, at x = y = ^ = 0. Note that the origin of the x-y-z 
coordinate system is arbitrarily selected, and other points may be used upon accounting for 
any translation. The sensor modules 120 - 126 receive the acoustic signals 165 from which 
they determine the acoustic range of the respective acoustic receiver 205 to the acoustic 
sources 110, 115. 

in a first embodiment the acoustic leceiver's position is determined analytically. In 
determining the acoustic ranges, the one way travel time ti for an acoustic signal 165 emitted 
ftom one of the high fi:equency, acoustic sources 1 10, 1 15 at position (xs, yj, Zs) to each sensor 
module 120 - 126 at position (x/, y,-, zd is given by: 

{^i-^sTHVi- y, )' + - 2 J = ^otf (1 ) 

Eq. (1) can be repeated for different acoustic sources, and, thus, for each of the 
acoustic sources 1 10, 115. Note that Eq. (1) constrains the position of the sensor modules 120 
- 126 to within a sphere centered on (x/, yu zd with radius a v^. The distance dr between the 
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sensor modules 120 - 126 is fixed (e.g.^ 12.5 m» in one particular embodiment) and can be 
approximated, assuming lixnited curvature of the cable 130 between sensors, by: 

Um - + [ym - y/f + - ^/P = dr^ (2) 

The inline angle of the sensor housing 215, or angular orientation of the acoustic 
receiver 205, with die horizontal plane is a derivative of the position vector. This can be 
represented by: 

7 ^^frr^ ^ = tan^^, (3) 

Eq, (2) and Eq. (3) can then be combined to constrain the positions of the sensor 
modules 120 - 126 to a plane. 

Various other measurements can also contribute to the solution in addition to, or in 
lieu of, the angular orientation: the water depth of the sensor modules 120 - 126; the distance 
between the vessel 105 and the position of the closest sensor; a depth measurement derived 
£com the hydrostatic pressure; measurements taken with a second acoustic source 115; and 
magnetometers (not shovm) inside the seismic cable 100. Note that these are all non-acoustic 
constraints. Li embodiments employing these other non-acoustic constraints, the sensor 
modules 120 — 126 will contain appropriate sensors in addition to, or in lieu of, the angular 
orientation measurement devices 210. In one particular embodiment, the water depth of the 
position of the acoustic receiver 120 - 126 is determined. One implementation of this 
particular embodiment measures the depth with a pressure sensor. A second implementation, 
however, retrieves the information from a data archive depending on other sensed data. Thus, 
the ascertaitmient of die non-acoustic constraint may include, e.g.^ retrieving an archived 
water depth measurement for the acoustic receiver's position. 

Note that, when an acoustic source acoustic sources 110, 115) is located on the 
vessel 105, it will be inline with the cable direction, but the cross-line direction is poorly 
constrained. Only a second acoustic source offset from the inline direction and angular 
orientation measurement devices 210 inside the seismic cable 100 can resolve this cross-line 
issue. Thus, in practice, the invention employs at least the two acoustic sources. Alternative 
embodiments, however, may use three or more acoustic sourx^es to facilitate the position 
determination. Theoretically, however, a single source may be employed in some 
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embodiments to generate separate acoustic signals at two different signal source positions 
from which two different acoustic ranges may be obtained. 

Thus, the position of each sensor module 120 - 126 can be determined analytically 
from their angular orientations constrained by, e.g., acoustic travel times. More particularly, 
the positions can be computed as the intersection of two spheres and a plane. FIG. 4 
conceptually illustrates the analytical solution 400 which is the intersection 405 of two 
spheres 410, 41S and a plane 420. Referring to both FIG. 1 and FIG. 4, the travel times from _ 
tfie two acoustic sources 110, 115 to a given one of the sensor modules 120 - 126 define the 
two spheres 410, 415 and the plane 420 represents the dip measurement resulting from the 
angular orientation of the respective sensor module 120 - 126. Each of the sensor modules 
120 - 126 is separated by the distance dr on the seismic cable 100, as was mentioned above. 

Generically, there are two points of intersection 405, shown in FIG. 5. In practice, 
using a third acoustic source 165 (shown in FIG. 1), this can be reduced to one point of 
intersection 405. Thus, the position of each sensor module 120 - 126 can be uniquely 
determined. To analytically determine the position for any one sensor module 120 - 126, the 
position of one acoustic source 110, 115 is temporarily assigned to be the origin. Each of the 
spheres 410, 415 is then represented by: 

(x,x) = r^ (4) 

and: 

(x-"m,x-m) = s^ (5) 
where m is the location of the second acoustic source 110, 115 relative to die first 
acoustic source 1 10, 1 15 and x is the position of the sensor module 120 - 126. 

The equation for the plane 420 is: 

(X4i) = d (6) 

where n = (0,0,1). The intersection of the two spheres 410, 415 is a circle 500, shown 
in FIG. 5. This circle 500 is confined to a plane given by: 

(x,m) = ((m,m) + r^ -s.^)/2 (7) 
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The intersection of this plane with the *dip* plane 420 is a line 420 with a 
parameterization a+X(nxm). Define a to be a linear combination of n anjd m: 

a=xn+ym (8) 

where 

(n,n)r-(n,m)s 
(n,n)(m,m)-(n,m)^ 



and 



(n,n)s-(n,m)r 
(n,n)(m,m)-(n,m)^ 



Ttie points of intersection of the line 420 with the first sphere 410 ^ves us two values 

for X: 



\^ = ±V(r"-(a,a))/(v,v). (11) 
Tlie two possible positions of the hydrophone now are a + A^v, anda+AjV, 
represented by the points 405 in FIG, 5. 

Using the measurement from the third acoustic source 165, or another constraint, one 
of the locations intersections 405) can be selected that corresponds to the trae position of 
the sensor module 120 - 126. Having deteraiined the position of the first sensor module 120 
- 126, the location of the second sensor module 120 - 126 can be determined in like manner, 
and so on. Eventually, the position of each sensor module 120 - 126 can be determined. The 
manner in which this additional constraint is applied will depend, to sonae degree, on the 
constraint itself. For instance, in one embodiment, one intersection 405 can be eliminated as 
an improbable physical location because it is, e.g., above the surface 150 or ahead of the 
survey vessel 105. 

The determination of the acoustic receiver's position (at 330, in FIG. 3) by the 
analytical technique discussed immediately above is represented in BIG. 6. The analytical 
determination includes first calculating (at 610) the intersection (e,g., the intersection 405) of 
a first sphere (e.g., the sphere 410), a second sphere (e.g., the sphere 415), and a plane (e.g., 
the plane 420). The first sphere and the second sphere are defined by the positions of a first 
acoustic source (e.g., the acoustic source 110) and a second acoustic source (e.g., the acoustic 
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source 115), respectively, relative to the position of the respective sensor module 120 - 126. 
Next, the analytical deteimination (at 1010) selects (at 620) one point of the intersection to 
identify the position of the respective sensor module. The illustrated embodiment 
accomplishes this by applying another constraint, e.g., acoustic travel time from a third 
acoustic source (e.g., the acoustic source 123). 

A second embodiment determines the position of the acoustic receiver through a 
model-based approach using an inversion algorithm. Typically, there are errors in the data 
(both the travel time and the dip measurements). These errors are not taken into account in 
the analytical determination of the position of the seismic cable 100. These errors are 
expected to propagate, at least to some degree, because of the iterative nature of the illustrated 
embodiment. To determine the effect of these errors, this second embodiment applies a linear 
inversion method to correct for these errors. 

More particularly, Eq. (1) -Eq. (3) are linearised by perturbing the position vector, as 
in Eq, (12) - Eq. (14) below: 

x^=xl+dxi (12) 

z,=z,+dz, (14) 
Ihseiting Eq. (12) - Eq. (14) into Eq. (1), and dropping second order teems gives: 

In a similar way, tiie substitution of Eq. (12) - Eq. (14) into Eq. (2) gives: 

2[x, -x,«]cfa, +2[x,« +2[y, -y,„]flfy, +2[y,« - y,]^fy,+i +2[z, -z,^,]ife, +2[«,^, -2, 

DifEraentiation of Eq. (3) then gives Eq. (17): 
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(z -zY 

Adx, - Adx,^i + Bdy, - Bdy,^^ + Cdz^ - Cdz,« = tan V, -7 ^ (17) 



[(^m-^,)'+(ym-y/)'f 



(18) 



B = 



2{yM-y,){zM-z,y 

hM'^if'+iyM-ytff 



C = i -2(gm-gi) ^ 

[{xM-Xif+{ym-yif\ 



(19) 



(20) 



For two sensors and two acoustic sources these equations can be written as a linear 
matrix equation Eq. (21): 



2(pc,-x^,) 2(y,-yJ 2Cz,-zJ 

%x,-x^ 2(ji-y^ Z.z,-z^ 

0 0 0 

0 0 0 

2^-2^« 2y-2yM 2zi-2zi« 

A -A B 



0 0 0 

0 0 0 

2(^1-%) 2(3i+,-yrf) 2(2J«-^,) 

2(^+1-%) 2()^44-3'J 2(2i«-5i2) 

2^«-2*} 2xi„-2^ 2x;„-2x; 
-J? C -C 



-(^ -^.if -yJ -{z, -z,if 

^ -U -^f -{yi -yJ -^f 
"^^-{xm-xj-bm-yj-iztn'zj- 

-(^ -^^7 -b+i -yJ" -(^« -^f 



taif 



(21) 



When three acoustic sources are used two additional lines similar two the top two lines 
are added. 
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A solution to the non-linear equations Eq. (1) - Eq. (3) is found, in the illustrated 
embodiment, using the iterative non-linear inversion technique taught by Tarantola, A. and 
Valette, B., "Generalized Non-linear Liverse Proljlems Solved Using the Least Squares 
Criterion," Rev. Geophys- Space Physics, 20, 219-232 (1982). This technique applies Eq. (9): 



where: 

po ^ the initial position vector of the cable; 
Pk = the position after the iteration; 
^PoPo - * he model covariance matrix; 

5 the data covariance matrix; 
do = the data; 

g s the forward model operator; and 

Gk = the matrix containing the partial derivatives (see equation (23) of Tarantola and 
Valette (1982)). 

The covariance and resolution matrix can also be computed. The term do-g(Pk) in Eq, 
(22) is equal to the right-hand side of Eq. (21), Gk equal to the left-hand side matrix in Eq. (2), 
and Pk to the left-hand side vector. 

Note that, although the illustrated embodiment employs the iterative non-linear 
inversion technique of Tarantola and Valette, oth» suitable techniques known to the art may 
be employed. It^tive non-linear inversion techniques are well known to the art and many 
are readily and commercially available in off-the-shelf computational software applications. 
For instance, the MatLab™ software package typically includes suitable techniques. 
Alternative embodiments might also employ other techniques, such as least-squares fit. 

The determination of the acoustic receiver's position (at 330, in FIG. 3) by the model- 
based technique discussed immediately above is represented in FIG. 7. Tlie method 700 
begins by modeling (at 710) the acoustic receiver's position from historical positions 
associated with the acoustic receiver's position. Next, the method 700 applies (at 720) an 



P.+i=Po + C^-G/-(c^+G,. 
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inversion algorithm to constrain the modeled position with the non-acoustic constraint. The 
method 700 is then iterated as the position of the acoustic receiver changes over time. 

Many aspects of the present invention are, in the illustrated embodiments, 
implemented in software, although the invention is not so limited. In altemative 
embodiments, these aspects may be implemented in electronic hardware or some combination 
of hardware and software. Nevertheless, some portions of the detailed descriptions herein are 
presented in terms of a software implemented process involving symbolic representations of _ 
operations on data bits within a memory in a computing system or a computing device. These 
descriptions and representations are the means used by fliose in the art to most effectively 
convey the substance of their work to others skilled in the art. The process and operation 
require physical mianipulations of physical quantities. Usually, though not necessarily, these 
quantities take tiie form of electrical, magnetic, or optical signals capable of being stored, 
transferred, combined, compared, and otherwise manipulated. It has proven convenient at 
times, principally for reasons of common usage, to refer to these signals as bits, values, 
elements, symbols, characters, terms, numbers, or the like. 

It should be bome in mind, however, that all of these and similar terms are to be 
associated with the appropriate physical quantities and are merely convenient labels applied to 
these quantities. Unless specifically stated or otherwise as may be apparent, throu^out the 
present disclosure, these descriptions refer to tiie action and processes of an electronic device, 
that manipulates and transforms data represented as physical (electronic, magnetic, or optical) 
quantities within some electronic device's storage into other data similarly represented as 
physical quantities within flie storage, or in transmission or display devices. Bcemplary of the 
terms denoting such a description are, without limitation, the terms ••processing," 
"computing," "calculating," "determining," "displaying," and the like. 

Note also that the software implemented aspects of the invention are typically encoded 
on some form of program storage medium or implemented over some type of transmission 
medium. The program storage medium may be magnetic (e.g., a floppy disk or a hard drive) 
or optical (e.g., a compact disk read only memory, or "CD ROM"), and may be read only or 
random access. Similarly, the transmission medium may be twisted wire pairs, coaxial cable. 
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optical fiber, or some other suitable transmission medium known to the art. The invention is 
not limited by these aspects of any given implementation. 

The vessel 105 is therefore equipped with a rack-mounted computing apparatus 800, 
conceptually illustrated in FIG. 8, with which these aspects of the invention executed The 
computing apparatus 800 includes a processor 805 communicating with some storage 810 
over a bus system 815. The storage 810 may include a hard disk and/or RAM and/or 
removable storage such as a floppy magnetic disk 817 and an optical disk 820. The storage 
810 is encoded with a data structure 825 storing the data set acquired as discussed above, an 
operating system 830, some user intetface software 835, and an application 865. The user 
interface software 835, in conjunction with a display 840, implements a user int^ace 845, 
The user interface 845 may include peripheral I/O devices such as a keyboard 850, a mouse 
855, or a joystick 860. The processor 805 runs under the control of the operating system 830, 
which may be practically any operating system known to the art. 

The processor 805, under the control of the operating system 830, invokes the usct 
interface software 835 on startup so that the operator can control the computing apparatus 
800. The application 865 is invoked by the processor 805 under the control of the operating 
system 830 or by the user through the user interface 845. The application 865, when executed 
by the processor 805, determines the shape of the cable either analytically, as in FIG. 4 — 
FIG. 6 or by the iterative, model based approach as in FIG. 7, depending on the embodiment 
implemented. The application 865 may also perform other functions, such as displaying the 
determined position. 

The invention may find many applications in the context of seismic surveying* In the 
illustrated embodiment, for instance, the invention may be used to determine the position of 
each of the acoustic receivers 205 for each of the sensor modules 120 - 126 where they are 
deployed for a survey. As will be appreciated by those skilled in the art having the benefit of 
this disclosure, that knowledge is, in itself, useful in conducting a seismic survey. However, 
because multiple positions are known, and because they are constrained to the cable 130 — ^in 
the illustrated embodiment — the shape of the seismic cable 100 can be determined. 
Furthermore, if the positions are determined dynamically during deployment, the shape of the 



SUBSTITUTE SHEET (RULE 26) 



wo 2004/036254 PCT/GB2003/004476 

16 

seismic cable 100 can be used to steer it during the deployment to help position the seismic 
cable 100 where desired for the seismic smrvey. 

For instance, FIG, 9 illustrates a method 900 by which the present invention may be 
implemented for dynamically detennining the shape of a body, e.g., the seismic cable 100. Jn 
general terms, the sensor modules 120 - 126 take certain measurements from which the 
effects of the deviations during deployment can be determined. From the known deviations, 
the position of the sensor modules 120 - 126 can be determined. In tum, the shape of the 
seismic cable 100 can be determined from the known positions. More particularly, the 
method 900 begins by measuring the dynamic angular orientation of at least two points on the 
body (at 910). This may include, for example, measuring the inclination and roll at fliat point 
on the cable. Note that the phrase "dynamic angular orientation," as used herein, implies that 
the point is moving or is subject to movement The method 900 then determines the position 
of the two points from the measured angular orientations dynamically (Le., on-the-fly) (at 
920). Note here that the term "determination" encompasses a determination within some 
acceptable degree of error, as no such determination is without some degree of inaccuracy. 
The term "dynandcally" implies in real-time or in near real-time. 

In one embodiment, dynamically determining the position of the two points from the 
measured angular orientation (at 920, in FIG. 9) comprises a method illustrated in FIG. 10. 
The method 1000 comprises first analytically determining (at 1010) the respective positions 
of a plurality of points (e.g., the sensor modules 120 - 126) on a cable (e.g.^ the seismic cable 
100). In the context of the present invention, flie phrase "analytically determined" means to 
calculate from actual measurements, as opposed to projected or predicted positions. One 
technique is discussed above relative to FIG. 4 - FIG. 6. In the that embodiment, these 
measurements are of inclination, roU, acoustic travel time (or range), although some 
implementations further include, e.g., water depth and/or heading measurements. From these 
analytically determined positions, the method 1000 determines (at 1020) the shape of the 
cable. 

In another embodiment, as shown in FIG. 11, a model-based method 1100 is 
employed to determining the position of the two points from the measured angular orientation 



SUBSTITUTE SHEET (RULE 26) 



wo 2004/036254 PCT/GB2003/004476 

17 

(at 920, in FIG. 9). The method 1100 is one implementation of the model-based technique 
previously discussed relative to FIG. 7. The method 1100 begins by modeling the respective 
positions of a plurality of acoustic receivers as the seismic cable 100 is deployed (at 1110). 
Any conventional modeling package known to the art suitable for this purpose may be used. 
The method 1100 then applies an inversion algorithm (at 1120) to constrain the modeled 
positions with the respective measured dynamic angular orientation of the respective 
positions. The non-linear inversion technique discussed above is applied at predetermined 
tune intervals to include the uncertainties and update the cable position. The cable position 
(i.e., the positions of the sensor modules 120 - 126) from the previous inversion is used as a 
starting point from which to find the small changes to determine the new position. Such a 
process mininoizes the number of iterations in the inversion process and speeds the delivery of 
the solution. The fiill solution is nevertheless periodically fully determined (at 1130) as a 
check to ensure that systematic errors are not accumulating in the real-time solutions, and that 
the error function for the solution is indeed the global minimum. Hie non-linear inversion is 
then applied (at 1130) to this new analytically determined position. 

A number of simulations were performed to quantify the value of adding additional 
measurements to the inversion for cable position. Each simulation was performed using a 
seismic cable including ten sensor modules, dipping from the horizontal at an angle of 22.50**, 
and with a sensor module spacing of dr = 12.5 m. The position of the first sensor was (1 1.54, 
0, -4.78). In all the simulations, the assumed initial position of the cable was determined 
analytically as described above. This initial position is exact in the case of zero errors in the 
data. Errors in the data were accounted for by performing a non-linear inversion, as was 
discussed above. This non-linear inversion can be repeated until convergence is obtained. 
The new cable position estimate can be compared with the known exact final position of the 
cable to find the accuracy, or error, in the estimated position. These errors are presented for 
each simulation- Failure to include the uncertainties in the inversion process may lead to the 
introduction of systematic errors in the seismic cable's position estimate. The included 
uncertainties are listed in Table 1 below. These simulations show that including dip 
measurements can significantly reduce the error in position of the seismic cable. 

Table 1. Measurement Uncertainties In Simulations 
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Parameter 



Uncertainty 



Sensor dip 



0.5" 



Water velocity 



30m/s 



Distance between 
sensors 



0,05m 



In a first simulation, only one acoustic source and no dip measurements were used 
The source position was (-12.5, 0, 0). The system of linear equations is now underdeteraiined 
and regularized using the data and model covariance matrices. Plots of the cable position 
before and after inversion are shown in FIG. 12, which projects the cable position onto the jc- 
y plane (top) and x-z plane (bottom) for this simulation. The true cable position is represented 
by the curves 1200, 1205, the perturbed cable position by the curves 1210, 1215, and the 
position after inversion by the curves 1220, 1225. The errors {Ex, Ey, in the position of 
each sensor and the square sum {E^ = sqrt{El + JSj + E^)) of the linearised solution are 

shown in FIG. 13. On average, E^^rz = 14 m for one acoustic source witti no dip measuremratt. 
For the second simulation, three acoustic sources in a triangular shape were positioned at (- 
12.5, 0, 0), (37.5, 70, 0), (37.5, -50, 0). Plots of the cable position are shown in FIG. 14. The 
em>is (JExj Ey^ E^ in the position of each sensor are shown in FIG. 15. The additional two 
acoustic sources reduces (he error particular in the crossline direction to Ex^r=^ m - 10 m. For 
the third experiment, dip angle measurements along the cable and the three acoustic sources 
were included. Plots of the cable position are shown in BIG. 16. The errors (£;!:, £y, iSy in the 
position of each sensor are shown in FIG. 17. The dip measurements further reduce tiie error, 
particular in the z-coordinate, to £:^z= 2 m - 4 m. 

Note that, although the illustrated embodiment pertains to a seismic cable deployed in 
water, the invention is not so limited and may be employed in other contexts. The mediod of 
the invention can be extended to the deployment of any cable in any medium, or even in a 
vacuum. Note that there is no such thing as a true vacuum, and even environments considered 
a "vacuum" contain some form of medium. Some of these variations may affect the 
availability, desirability, or selection of constraints that may be used in determining the 
positions from the angular measurements. Similarly, the invention may be applied to 
stationary cables, a cable previously deployed and at rest on the sea floor. In these 
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embodiments, the angular orientation measurement and the detemiination need not be made 
dynamically and the determination need not be made on the fly. These and other variations 
on the embodiments disclosed herein will become apparent to those skilled in the art having 
the benefit of this disclosure. 

Thus, the invention uses, in the illustrated embodiment, combination of angular 
measurements made along the length of a cable with other constraints such as, acoustic arrival 
times between acoustic sources and acoustic sensors also mounted on the cable, to determine 
the cable's shape in real-time. The angular measurements provide additional constraints that 
can improve the accuracy of a solution based on acoustic measurements only. It may also be 
possible to reduce the complexity of the acoustic system through the use of additional 
information afforded by measurements of angle. In the illustrated embodiment, the technique 
disclosed herdn is implemented such that cable position can be determined within 2-3 
seconds of measurements being taken, in pseudo real-time. This permits one to effect control 
of the seismic cable's shape in the catenary as it is being deployed. This real-time, or pseudo 
real-time, cable position determination (including uncertainties) lends itself to an iterative 
process. 

The various embodiments and aspects of the invention disclosed herein consequentiy 
include a method and an apparatus for determining a position of an acoustic receiver. The 
apparatus includes at least one acoustic source; an acoustic receiver, and a computing system. 
TTie acoustic receiver is capable of receiving a plurality of acoustic signals transmitted by the 
at least one acoustic source from at least two signal source positions. The computing system 
is programmed to determine a position of the acoustic receiver ftom the acoustic ranges 
between the at least two signal source positions and the acoustic receiver and a non-acoustic 
constraint The method includes determining a first and a second acoustic range from a first 
signal source position and a second signal source position, respectively, to the acoustic 
receiver, ascertaining a non-acoustic constraint on the acoustic receiver's position; and 
detennining the acoustic receiver's position from the first and second acoustic ranges and the 
non-acoustic constraint. 
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This concludes the detailed description. The particular embodiments disclosed above 
are illustrative only, as the invention may be modified and practiced in different but 
equivalent manners apparent to those skilled in the art having the benefit of the teachings 
herein. Furthermore, no limitations are intended to the details of constraction or design herein 
shown, other than as described in the claims below. It is therefore evident that the particular 
embodiments disclosed above may be altered or modified and all such variations are 
considered within the invention as claimed below. Accordingly, the protection sought herein 
is as set forth in the following claims. 
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